Pseudogap in cuprate and organic superconductors 
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We study the pseudogap present in cuprate and organic superconductors. We use the dynamical 
cluster approximation (DCA), treating a cluster embedded in a bath. As the Coulomb interaction 
is increased, cluster-bath Kondo states are destroyed and bound cluster states formed. We show 
that this leads to a pseudogap. Due to weaker coupling to the bath for the anti-nodal point, this 
happens first for this point, explaining the k-dependence of the pseudogap. The pseudogap can be 
understood in terms of preformed d-wave pairs, but it does not prove the existence of such pairs. 

PACS numbers: 71.10.-w; 71.27.+a; 71.10.Fd 



Experiments show a pseudogap at the Fermi energy 
in the normal phase in both cuprate[l| and organicQ su- 
perconductors under certain circumstances. For cupratcs 
the pseudogap forms for k = (ir, 0) while for k = 
(tt/2, tt/2) the spectrum has a peakQ. It is crucial for 
the understanding of these systems to trace the origin of 
the pseudogap. 

Calculations using embedded cluster methodsQ, e.g., 
the dynamical cluster approximation (DCA) reproduce 
a k-dependent pseudogap for the Hubbard model[6- 13 1. 
Ferrero et a/.[14| discussed small embedded clusters in 
terms of a transition between a state where the cluster 
orbitals form Kondo-states with the bath and a state 
where the cluster forms a bound state and a pseudo- 
gap. However, the strong k-dependence of the pseu- 
dogap in cuprates was not discussed. More recently, 
the pseudogap has been interpreted as a momentum- 
selective metal- insulator transition ll[ . There has been 
much work relating the pseudogap to preformed super- 
conducting pairs [l5(, which have have not reached phase 
coherence and superconductivity at the temperature T 
studied. On the other hand it has been argued that the 
pseudogap and superconductivity phases competeflij. 

The DCA treats a cluster of N c atoms embedded in 
a bath. Guided by the smallest (N c = 4) cluster giving 
a pseudogap, we construct a very simple two-site two- 
orbital model. For a small Coulomb interaction U, the 
cluster orbitals form Kondo states with the bath. As 
U is increased, the Kondo states are destroyed and a 
bound state is formed on the cluster. We show that 
if this state is nondegenerate, this leads to a pseudo- 
gap. By comparing correlation functions, we find that a 
N c = 8 DCA calculation behaves in a similar way. Due 
to the weak dispersion around k = (71", 0), the coupling 
of this cluster k-vector to the bath is much weaker than 
for k = (tt/2, tt/2). Then the Kondo state for (tt/2, tt/2) 
is destroyed for a larger U, giving the k-dependent pseu- 
dogap. We show why the pseudogap is lost for a frus- 
trated electron-doped cuprate. We find that the pseu- 
dogap can be interpreted in terms of d-wave supercon- 



ducting pairs. The pseudogap hinders the propagation 
of pairs and hurts superconductivity. 
Below we study the Hubbard model 
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where c\ a creates an electron on site i with spin a , ni a = 
c\ a c ia , tij = t if i and j are nearest neighbors, t' if they 
are second nearest neighbors and zero otherwise. Here 
we use t = —0.4 eV and t 1 = or — 0.3i. U is a Coulomb 
integral and \x the chemical potential. 

We solve this model using DCA, and discuss the cluster 
in terms of its N c K-states. With N c — 8 there are both 
nodal [K = (±tt/2, ±tt/2)] and antinodal [(tt, 0), (0, tt)] 
states. For isolated clusters, the nodal and antinodal 
spectra are identical for N c — 8 and 16. This is due 
to a symmetry for these values of iV c Q. Since DCA cal- 
culations with N c = 8 and 16 give a pseudogap for the 
antinodal point, this must then be due to the coupling 
to the bath. To test this we have switched the baths. 
The pseudogap then indeed appears at the nodal point 
instead of the antinodal point. An important clue is that 
the coupling of the bath to the antinodal point is weaker 
by a factor of three to four, due to the weaker disper- 
sion of the band around this point. For a self-consistent 
DCA calculation, we find that the bath also tends to have 
a pseudogap. Is then the pseudogap in the spectrum 
mainly due to the pseudogap in the bath or to cluster 
properties? We find that performing only one iteration, 
with a fully metallic bath, the pseudogap is reproduced, 
although for a larger U value. Below we then focus on the 
first iteration. The key issue is then to understand the 
associated multi-orbital quantum impurity. We assign 
the pseudogap to the weakness of the antinodal coupling 
and internal electronic structure of the cluster. 

The smallest embedded cluster with a pseudogap is 
N c = 4. K = (0, 0) is mainly occupied and K = (tt, tt) 
mainly unoccupied. We can then focus on K = (tt, 0) 
and (0,7r), which mainly determine the dynamics. For 
simplicity, we first only couple each K-state to one bath 
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state. This leads to a two-site model with a two-fold 
orbital degeneracy. This model has the essential features 
of a pseudogap and we can study the origin in detail. 
We introduce a hopping integral V between each cluster 
orbital, labeled lc or 2c, and its bath site orbital, labeled 
lb or 2b. On the cluster site c there is a direct Coulomb 
integral U KX = U — AU between spin up and down states 
of the same orbital, U xy = U + AU between two different 
orbitals and an exchange integral K. The lowest S~ = 
states on the isolated c site are 

I IT) = ^(4c t c U T 4 t c U); Ei* = ^xx T K (2) 
l 2 T> = -^{c\ c Aci T 4t c ici); E 2T = U xy T K, 

where cj c . creates a spin up electron in level i on site c 
and E izfl are the energies of the states. The state |2— ) 
is a triplet and the other states are singlets. To simulate 
the lowest three states of an isolated cluster with N c = 
4, we identify lc with (n,0) and 2c with (0,7r) and put 
AU = 0.03U and K = 0.117. The lowest state in the 
N c = 4 cluster and the model cluster are both singlets. 
Uxx < U xy may seem surprising, but is related to the 
influence of the (0, 0) and (71", tt) orbitals, neglected in the 
model cluster. Alternatively, we can choose f7 xx > f7 xy . 
Then the lowest model cluster state is a triplet and, we 
will see, the physics is completely changed. 

For the two-level, two-orbital model we have calculated 
the correlation function Sidb = (Si c • Sib), between the 
spins in orbitals lc and lb, on the cluster and in the bath. 
We have —3/4 < 5i c ib- This is shown in Fig.Q]as a func- 
tion of U (full curve) . For moderate values of U there is a 
strong negative correlation, which increases with U . The 
correlation function C\ c ^i c ^ = (ni c -|-ni c .|,) — 1/4 is also 
shown. It becomes negative as U is increased, implying 
a reduced double occupancy of orbital lc and a spin 1/2 
state starts to form on lc. Both these results are con- 
sistent with levels lc and lb forming a Kondo-like state. 
Further increase of U , however, leads to a rapid reduc- 
tion of |Sicib| and Ci c fi c ± grows positive. The Kondo 
singlets are then broken, and instead a bound state is 
formed on the cluster. Increasing U reduces the Kondo 
energy, while the gain from correlating the electrons on 
the cluster increases. In the lowest cluster state Cictiai, 
takes the value 0.25, which is approached in Fig. [TJ 

Fig.[T]also shows weights of the spectral function on the 
photocmission side for the two-orbital, two-site model. 
For U = all weight is in the leading peak. As U is in- 
creased, the weight of the Hubbard peak increases, cor- 
responding to final state with primarily one electron less 
on the cluster site. At the same time three peaks close 
to the leading peak get more weight, summed up in the 
curve "Satellites" . These peaks correspond to final states 
with substantial weights in excited neutral cluster states, 
while the leading peak corresponds to a final state with a 
large weight of the lowest neutral state. For t//|V| ~ 11 
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FIG. 1. Weight of the correlation functions 5i c ib and n lc ^ni a _i 
as well as spectral weights for the leading peak, satellites form- 
ing the peaks surrounding the the pseudogap and the Hub- 
bard peak for the two-level two-orbital model. 



the satellites acquire more weight than the leading peak 
and a pseudogap starts to develop. This happens when 
the Kondo-like effect starts to be lost and a bound state 
on the cluster starts to form. 

To address the pseudogap, we study a general clus- 
ter with U » \V\. Assume an integer number of elec- 
trons per cluster, no, and project out the piece, |$n )j 
of the ground-state, |$), that has exactly no electrons. 
We study photocmission processes, removing an electron 
with the quantum number v and energy e. We split |$ no ) 
as 

l*»o) = |0> + |1>, (3) 

where c„|0) = and cjcjl) = |1). Let H a = 
^2 VE {Vvecicvs T h.c.) give the hopping into the cluster. 
We then write the part of the ground-state corresponding 
to no + 1 electrons as 

v' e 

where we have approximated the energy difference be- 
tween cluster states with one extra electron (hole) and 
the lowest neutral state as AE+ (AE-). If the lowest 
state of the isolated cluster is nondegenerate, this is (non- 
degenerate) perturbation theory. However, if the lowest 
state is degenerate, e.g., for the Kondo problem, |0) and 
|1) contain much more information than can be obtained 
from perturbation theory, and the treatment goes beyond 
perturbation theory. For an infinite system, we construct 
a state \vef) with a hole in the bath at the Fermi energy 
£f- We form c„ £F |$ no ) and allow hopping to states with 
no — 1 electrons on the cluster. The corresponding am- 
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FIG. 2. Initial and final states for a process removing an elec- 
tron with quantum number v. There is a negative interference 
between two channels for reaching the final state. 



plitude is then 
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If the isolated cluster ground-state is nondegenerate and 
AE± are large, the number of holes in the bath for the 
states |0) and |1) is small. Then the third and fourth 
terms are also small, while the expectation values in the 
first and second terms are large and tend to cancel. If the 
orbital v is half- filled the cancellation becomes perfect. 
This cancellation is illustrated schematically in Fig. [2] 
For, e.g., U = 1.0 and V = —0.02 in the model above, 
the four terms are -0.0117, 0.0121, 0.0020, 0.0020. This 
illustrates the strong cancellation between the first two 
terms and the small magnitude of the next two terms. 

If the lowest cluster state is degenerate, the situation is 
quite different. The system may then form a Kondo-likc 
state even for large U. The first term in the amplitude 
above is then smaller and docs not cancel the second 
term. Furthermore, the "spin-flip" terms three and four 
can be substantial. As a result, a pseudogap may not 
develop. In the model above, a negative AU gives a 
triplet ground-state for the model cluster. In the limit of 
a large U, there are then triplet states on both the c and 
b site, coupling to a total singlet. The leading peak has 
more weight than each of the neighboring satellite for the 
U range shown in Fig. [TJ and a pseudogap does not form. 
For an infinite system, these results can also be discussed 
in terms of the phase shift at the Fermi energy 14 1 . 

Wc now study N c = 8 at half-filling. Fig. [3] 
shows correlation functions of the type C , ( Wi0 )t,(7r.o)4. = 
(«(7r,0)t n (7r,0H) - («(ir,o)f ) ( n (^m) ■ This mnction behaves 
in a very similar way as Ci c -|-ic4 in Fig. [T] As U is in- 
creased it first turns negative, indicating the formation 
of a Kondo-like state. Then it turns positive, indicat- 
ing the formation of a bound state on the cluster in the 
(it, 0), (0, 7r)-space similar as for the two-site two-orbital 
model. This is also supported by the strong negative 
correlation between (n, 0) and (0,7r) as is also found for 
|1— ) in Eq. At the same time a pseudogap develops 
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FIG. 3. Correlation functions as a function of U for an em- 
bedded eight site cluster in the first iteration. The parameters 
are t = -0.4 eV and T = 290 K. 



TABLE I. Bath-cluster spin correlations and hybridization 
energies from ED calculations. The parameters are t — —0.4 
eV, V(ir,0) = 0.08 eV and V(ir/2,ir/2) = 0.21 eV. 









(S z s z ) (V) 


u 




n 


(M) (M) 


2.5 


1.25 


1.0 


-.0020 -.0051 -.021 -.097 


3.2 


1.6 


1.0 


-.0012 -.0030 -.014 -.066 



in the (?r,0) spectrum, as in Fig.ffl C( w / 2 , 7r /2)t,(7r/2, 7 r/2)^ 
behaves in a similar way, but it stays Kondo-likc up to 
larger values of U before the (±7r/2, ±7r/2)-spacc is also 
used to form a bound state on the cluster. Therefore the 
(±7r/2, ±7r/2) pseudogap forms for larger U. The reason 
is that the coupling to the bath is a factor of three to 
four stronger for (tt/2, 7t/2) and it is favorable to keep 
the Kondo-likc state up to a larger U. 

To further support this, we have performed exact di- 
agonalization (ED) calculations for N c = 8 and one bath 
level at the chemical potential with the coupling V(K) 
to the cluster K-state. Wc study the correlation (S z s z ) 
between spins on the cluster and the bath and the hop- 
ping energy (V) between the cluster and the bath for a 
certain K. Table [J shows that both quantities are much 
larger for (tt/2, tt/2) and that both decrease rapidly with 
U. 

The energy gain when all electrons of the isolated clus- 
ter correlate is much larger than the sum of the en- 
ergies gained when correlation is only allowed in the 
(■7T, 0), (0,tt)- or (±7r/2, ±7r/2)-spaces. This favors a si- 
multaneous switch of both spaces from the Kondo-likc 
states to a correlated cluster state. However, the dif- 
ference in the coupling to the bath is so large that the 
switch, nevertheless, happens for different U in Fig. [3] 
This can be different for the frustrated case. 

We now consider the frustrated case, t' — — 0.M, for 
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TABLE II. Bath-cluster spin correlations and hybridization 
energies from ED calculations on N c = 8 frustrated clusters 
with U = 2.5 eV, t = -0.4, t' = -0.3t, V(tt,0) = 0.08 eV and 
V{tt/2,-k/2) = 0.21 eV. We also show J(K) [Eq. ©]. 
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different fillings and U = 8\t\ in the first iteration. For 
moderate hole doping the N c = 8 calculation shows a 
clear pseudogap. On the other hand, for electron dop- 
ing we find no pseudogap. We first study this using ED 
as above. The coupling to the bath is kept fixed, and 
n is varied by varying fx. Table [IT] shows how both the 
cluster-bath spin correlation and hopping energy increase 
for (tt, 0) as we go from hole- to electron-doping, while 
the change is small for (ir/2,ir /2). To understand this, 
we study the J(K) obtained from a Schrieffer-Wolf trans- 
formation 



J(K) = |V(K)|-(^ 



1 



(K) + U-n (J, — e(K) 



(5) 



Due to the shift of /i with n, J(K) increases for (w, 0), 
while the levels for (tt/2,tt/2) are located in such a way 
that the changes of J(K) are relatively small. 

In the ED calculation we assumed a doping in- 
dependent coupling. The coupling is related to Im 
l/Go(K,w n ), where Go is the Green's function of the 
bath. For t 1 = — QM and moderate hole doping, both 
Re Go[(tt, 0), uj n ] and Im Go[(7T, 0), u n ] are large for small 
imaginary frequencies w n . Then Im l/Go[(7r, 0), w n ] is 
small and the coupling relatively weak. As n is increased, 
/i shifts so that for the electron doped system both the 
real and imaginary parts of Go are smaller and the cou- 
pling larger. This increase and the increase of J(K) with 
n favor a Kondo effect for (tt, 0) for the electron-doped 
system. The Kondo effect then tends to be lost simulta- 
neously in the (71", 0)- and (±7r/2, ±7r/2)-spaces. This re- 
sult depends crucially on the sign of t'/t. For the organic 
superconductors, with a triangular lattice and at half- 
filling, wc find that the pseudogap goes away at t' ~ 0.6t, 
in agreement with experiment [2j, since J(ir, 0) has a min- 
imum at t' = for half-filling. An additional effect is 
that when the system is doped, the ground-state has a 
substantial mixture of a cluster state with one electron 
more or less. This state is degenerate, which is unfavor- 
able for a pseudogap and in agreement with the lack of 
a pseudogap for large doping. 

To understand the character of the pseudogap, we have 
calculated Pd = ^•(AjAD/iVc, where Aj is the operator 



TABLE III. Energy E, spin S. degeneracy Deg. and (/-wave 
superconductivity correlation Pd for low-lying states with 
S z — for N c = 8. Deg. refers to the degeneracy of S z = 
states. The parameters are U = 3, t = —0.324 eV and t' = 0. 
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for d-wave superconductivity. A factorized term was sub- 
tracted so that Pd — for U — 0. For the ground-state 
of the isolated cluster Pd is rather large (see Table llllj) . 
This state dominates the ground-state of the embedded 
cluster for relatively large U. The peaks surrounding 
the pseudogap correspond primarily to final states with 
a large weight of excited neutral states with S = or 1, 
for which Pd is substantially smaller. We can then think 
of the pseudogap as corresponding to the break up of a 
d-wavc singlet pair. 

This relation of the pseudogap to d-wave pairs does 
not contradict the observation that the pseudogap and 
superconductivity competeflij]. The occurrence of super- 
conductivity is determined by the Bethe-Salpeter equa- 
tion, and depends on the matrix (1 — x°r c ) becoming 
singular Here r c is the two-particle irreducible vertex, 
representing the interaction, and \° consists of products 
of two dressed Green's functions, describing propagation 
of pairs. x° 1S closely related to the density of states at 
the Fermi energy. As U is increased, T c grows and x° is 
reduced. This reduction becomes very rapid as a pseudo- 
gap forms 17J . This hurts the onset of superconductivity. 
A similar competition between \° an d r c was observed 
for alkali-doped fullerides(l~8j]. 

From the discussion below Eq. ([4]) we may expect to 
see a pseudogap if the lowest state of the cluster is non- 
degenerate, independently of Pd- To test this, we have 
applied an antiferromagnetic (AF) field J2i ^( n it — n i\)i 
where Si = 5 on one sublattice and —S on the other. This 
term is only applied to the cluster and not to the bath. 
One iteration is performed. This term favors an AF state 
on the cluster and reduces Pd- Fig. |4] shows that there 
is, nevertheless, a pseudogap for K = (tt, 0). 

In conclusion wc have provided and analyzed a very 
simple model, showing pseudogap formation. As U is in- 
creased, the system switches from a case where the clus- 
ter orbitals form a Kondo state with the bath to a case 
where a nondegenerate bound state is formed on the clus- 
ter. This leads to a pseudogap. We showed that this anal- 
ysis also applies to larger clusters. It is crucial that the 
coupling to the bath is much weaker for K = (tt, 0) than 
for (71-/2,71-/2), which leads to the strong K-dependence. 
We have found that the pseudogap is lost for the electron 
doped system, due to an increase in the effective J and 
the cluster-bath coupling for K = (71, 0). The pseudogap 
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FIG. 4. Spectrum averaged over k around K = (n, 0) and 
(it/2, 7r/2) when an AF field S is applied to the cluster. 



can be interpreted in terms of d-wave superconducting 
pairs. It strongly reduces x > describing how the propa- 
gation of pairs is hindered, leading to a competion with 
superconductivity. 
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